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Abstract 

We have established semiclassical kinetic equations for various spin corre- 
lated pumping phenomena incorporating adiabatic spin rotation in wave func- 
tions. We employ this technique to study topological pumps and illustrate 
spin pumping in a few models where various spin configurations or topolog- 
ical motors drive adiabatic pumps. In the Rashba model we find that a 
topological spin pump is driven by a meron with positive one half skyrmion 
charge, the size of which can be controlled by external applied gates or Zee- 
man fields. In the Dresselhaus model on the other hand, electron spins are 
pumped out by a negative meron. We examine the effects of Zeeman fields 
on topological spin pumping and responses of Fermi seas in various topo- 
logical pumps. The phenomena of topological pumping are attributed to 
the beam splitting of electrons in the presence of spin rotation, or topolog- 
ical Stern- Gerlach splitting and occur in a transverse direction along which 
charge pumping currents might either vanish or are negligible. The transport 
equations established here might also be applied to the studies of anomalous 
Hall effect and spin Hall effect as demonstrated in one of the appendices. All 
results are obtained in an adiabatic expansion where the adiabaticity condi- 
tions in Eqs. (36), (42), (99), (100) are satisfied. 
PACS number: 03.65. Vf, 73.43.-f 



I. INTRODUCTION 



Adiabatic transport of electrons or quantum pumping which is nearly reversible has been 
one of very promising means to manipulate coherent wave packets in the extreme quantum 
limit. In the presence of periodic adiabatic perturbations, a net charge can be transferred 
across a quantum structure during each period which is independent of frequencies of exter- 
nal perturbations and which represents a DC current induced by adiabatic perturbations. 
This phenomenon was first observed in an early work on transport of edge electrons in 
quantum Hall states [1]. But general solutions to this problem were provided in Ref. [2] 
where conditions of quantized charge transport were established. The robustness of quan- 
tized transport with respect to disorder potentials and applications to quantum Hall effects 
were later on studied in a series of works [3]. 

The absence of dissipation during the adiabatic process is evident if perturbations are 
applied to a closed quantum structure with nontrivial topology such as a mesoscopic ring 
or torus. If one further assumes that the electron spectrum is discrete and the external 
frequency is incommensurate with energy gaps in the spectrum, then no resonance absorption 
can occur and quantum states evolve via unitary transformation which conserves the entropy. 
The absence of entropy production in the adiabatic process is therefore a natural consequency 
of a pure state evolution which has been known for a while. 

For a quantum structure with continuous spectra either because of contact with leads in 
a mesoscopic limit or more generally because of thermal broadening, it is more convenient 
to introduce one-particle density matrices to describe the evolution of quantum systems. 
The adiabaticity can be achieved when the frequency of applied perturbations is lower than 
various relaxation rates characterizing the dynamics of one-particle density matrix. The 
issue of entropy production in this case however hasn't been fully addressed and is less well 
understood. Nevertheless, it is widely appreciated that dissipation involved in adiabatic 
transport in this limit should also be smaller than that due to a transport current with 
biased voltages applied. 
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Given the obvious advantage of the adiabatic charge transport, in this article I intend to 
generalize the idea to adiabatic spin transport. I will address the issue of spin pumping in 
both limits which can be easily achieved in laboratories and limits which are theoretically 
exciting but might not be as easy to be realized in solid state structures. Particularly, I will 
propose a novel spin pumping mechanism which is based on the topological beam splitting of 
electrons instead of the usual Zeeman splitting. Special classes of models are introduced to 
facilitate discussions on topological spin pumping. At the end, I will compare the efficiency 
of different spin pumping schemes. 

One obvious mean to pump spin out of the system is to adiabatically transfer polarized 
electrons in quantum structures. During the adiabatic transport, the currents carried by 
spin-up and spin-down electrons have an asymmetric part and therefore electrons pumped 
out of the structures also carry net spins. This standard scheme is reviewed in section III. 

In section IV, I discuss the phenomena of topological beam splitting in details. Par- 
ticularly, I demonstrate that spin rotation in either real-space (X-space) or in Fermi seas 
leads to transverse motion of electrons. In section V and VI, I investigate topological spin 
pumping due to spin rotation either in the X-space or in Fermi seas. In both cases, spin-up 
and spin-down electrons though both are electrically negatively charged, carry certain topo- 
logical charges with opposite signs. As consequencies, spin-up and spin-down electrons can 
be split because of opposite topological transverse forces, an analogy of splitting of electrons 
and positrons in an orbital magnetic field. I would like to refer this kind of beam splitting 
as topological Stern-Gerlach splitting ( TSGS) to contrast the usual Stern-Gerlach splitting 
of spin-up and -down particles in atomic physics. In this article I focus on the origin of this 
phenomenon and basics features. I plan to present a practical design of a topological spin 
pump in a subsequent paper. 

Finally, in connection with topological spin pumping to be discussed in the article, it 
is worth mentioning a few recent works on anomalous Hall effects and spin injection where 
the topology of Hilbert spaces plays a paramount role. In Ref. [7], the authors pointed 
out that interactions between conduction electrons and background skyrmion configurations 
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activated in magnetite might be responsible for the sign and temperature dependence of 
anomalous Hall effects observed in Ref. [8]. The authors of Ref. [9] meanwhile argued that 
k-space Chern-number densities also modify the equation of motion for electrons. The 
corresponding contributions to the anomalous Hall effect in ferromagnetic semiconductors 
were further studied in Ref. [10]. Following these works, it is now believed that the anomalous 
Hall effect can be an intrinsic phenomenon. It is indeed likely to occur when skew-scattering 
from impurity atoms is absent as proposed by Karplus and Luttinger a while ago [11]. In 
Ref. [12] the authors have considered intrinsic spin Hall currents in semiconductors; many 
interesting features have been found. Related discussions can be also found in Ref. [13]. 

Independently in a series of illuminating works [14-16] the authors studied spin injec- 
tion in semiconductors characterized by the Luttinger Hamiltonian. They have found that 
singular topological structures in the k-space as well have fascinating effects on accelerated 
electrons and lead to important consequencies on spin injection. In Ref. [15], the authors 
further pointed out possible connection between transverse spin Hall currents and supercur- 
rent in superconductors. The issue of dissipation however is still under debate and remains 
to be fully understood. 

II. KINETIC EQUATIONS FOR ONE-PARTICLE DENSITY MATRIX 

Consider the one-particle density matrix p Q , /3 (x / , x; t', t). Subsripts a, (5 — ± are intro- 
duced as spin indices; later in this article I also introduce 77, £ = 1,2, ..N as indices in an 
iV-dimensional parameter space; /i, v = x,y,z as indices in the real or momentum spaces. 
The evolution of one-particle density matrix is determined by the following equation 



[^ + ^]Pa/3(x',x;t / ,t) = H a /3'(x,i— ;t)p w (x',x;t / ,t) - p Q/3 /(x / , x; t', t)Hp'p(x', i-^j ; if). (1) 

To study the transport in a semiclassical limit, one introduces 

x + x' 



2 ' 



X = x' - x; 



T= t -±t,r = t>-t. (2) 



Furthermore, one defines a generalized semiclassical density matrix 

p a/3 (k,r;w,T) = 1 1 dXdrexp(ik.X-iwr)p^(r + y, r- *;T+^,T- T -) (3) 

V is the volume of systems. 

In a semiclassical approximation, one obtains the equation of motion for the one-particle 
density matrix, 

of + [ dk d~ r m dk + of aJ ] pMK r; w> T) 

= -// a ^(k,r;T)(g)p w (k,r;a;,r)+X c -V/3(k,r;a;,r). (4) 
Here 

Aap^Bpy = AapBpy — BapAp^. (5) 

and X C L is a collision integral operator for elastic (nonmagnetic) scattering processes, k 
and r in Eq.4 are variables instead of operators. The gradient expansion which is valid as 
far as the transport occurs at a scale much larger than the fermi wave length is sufficient 
for the study of semiclassical phenomena. In all models employed in this article I find the 
commutator H®p in Eq.4 vanishes in the semiclassical approximation. 

The charge current J and spin current J 2 with spin along e z direction are 
rTA _ fdu d 3 k dH a p(k,r) 

J ( r, 



' T) = i^JW dk ^(k,r;o,,T), 



2tt (2tt 

-tzi r r [ \ fducPk dH a ,p(k, r) 

) = J 2^(2^p z ' Gaa ' dk ^a(k,r;w,T). (6) 

To facilitate discussions on the adiabatic transport, one further separates the one-particle 
density matrix into symmetric and asymmetric parts (pf^), 

p a p(k, r; u, T) = pf /3 (k, r; u, T) + p^(k, r; u, T) 

l aL p%(k, r; T) = 0,I aL pUk, r; T) = lp^(k, r; T). (7) 



And 



J rfO(k)p^(k,r;cu,T) = J dil (k) H (k) pf ^ (k, r ; a; , T) = 0, (8) 

f2(k) is introduced as a unit vector along the direction of k in Eq.8 and in the following 
sections. 

In the relaxation approximation employed here, elastic nonmagnetic impurity scattering 
only leads to momentum relaxation because the collision integrals are SU (2) singlet operators 
and act trivially on the density matrix p a p. This however doesn't generally imply that 
impurity scattering combined with the spin-orbit coupling which I am going to discuss 
should not cause transitions between different spin states. Nevertheless in the adiabatic 
approximation employed in this article, in a special basis these transitions are negligible 
(see discussions about SU(2) gauge fields and adiabaticity conditions in section V and VI). 
As far as the adiabaticity conditions are satisfied, the collision integral can be treated in 
the usual Born approximation even in the presence of spin-orbit coupling. Please see more 
specific discussions about the adiabaticity at the beginning of section VI, and discussions 
after Eqs.72,100. 

Since I am interested in the transport phenomena at distance much longer than the mean 
free path l or at frequencies much lower than the inverse of mean free time r , i.e., 

L > l = t v f ,T > r , (9) 

I adopt the standard diffusion approximation. Furthermore for the study of adiabatic charge 
and spin pumping phenomena it is sufficient to keep the first order term in an adiabatic 
expansion. Taking into account the definition of symmetric and antisymmetric components, 
one obtains p A and p s as 

p% = Wt n„(k)[- J- V + m '-"'^- T 1 g, ^-K,„(k, r; u,, T) (10) 



Vk = |k|/m is the velocity and Dk = vlr /d is a diffusion constant; al = 2, 3 are the 
dimensions of the Fermi seas which interest us. p° al3 is the equilibrium one-particle density 
matrix. 
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The charge pumping current and spin pumping current with spin pointing at the z- 
direction can then be expressed as 

T2 f duo d 3 k r <9 . dH a > (k, r; T) 9 . <? ^ x 

J - = / a? (2^ D ' e ' ' ^' ["ft"*" + % S^ (k - r ^ T) ' 

(12) 

Here and in the rest of the article I set h — e — 1. This set of equation will be used to study 
various spin pumping phenomena. 

III. CHARGE AND SPIN PUMPING OF POLARIZED ELECTRONS 

I first apply the kinetic equations to study adiabatic charge transport of polarized elec- 
trons. The external perturbations are represented by N external a. c. gates with same period 
T ; i.e. 

V ext (r,T) = Y / 9v(T)V ri (r-r v ), 

v 

g v (T)=g v (T + T ),r ] = 1,2,. ..N. (13) 

For pumping phenomena, the boundary conditions at r = ±L/2e fM (/i = x,y) are chosen as 

P%(K "§e M ; u, T) = p%{k, u, T). (14) 

Eq.14 is valid when a) the sample has a closed geometry along e M or b)more practically leads 
at boundaries are ideal and are maintained in a thermal equilibrium, which also corresponds 
to a current biased situation. 

As noticed in a previous work [17], at time T the one-particle density matrix in the pres- 
ence of adiabatic perturbation only depends on the potentials at that moment. Particularly, 
the density matrix is a function of g n {T), rj — 1, 2, ..N and their time derivatives g n iT); and 
it has this local time dependence as a result of adiabaticity. The charge transport per period 
T therefore has the following appealing general structure (77, 77', ^ = 1,2, ...TV) [17], 



Qu = eJd 9ri A dg^,n r , = ( — — - __)Tr{— p(k, r; {^}, (15) 

Here I introduce Qa as the charge transport along the ej-direction. dg v A dg% is a skew 
symmetric wedge product, i.e., dg v A = — dg^ A dgv,. The trace Tr is carried over the 
momentum, real-space and spin space. In Eq.15, the charge transport has been expressed 
explicitly in terms of the adiabatic curvature n v ^; the form of the curvature is uniquely defined 
by the local time-dependence of the one-particle density matrix p a /j(k, r; {g v >(T)}, {g^iT)}). 

In the following I am going to evaluate the one-particle density matrix p aj g(k, r; T) and 
therefore the curvature explicitly including the spin polarization. As indicated in Eq.ll, 
the antisymmetric part of the density matrix can be expressed in terms of the symmetric 
part. And the symmetric part of the density matrix receives a nonadiabatic correction 
following the second line in Eq.ll; the solution is 

p%(k, r; u; , T) = p%(k; u) + ^-Mf^k, r;u,T), 

M3(k,r; W ,T) = J dv'G(v, v')g v ^-H a ^ r'; {g rf (T)}) ® U) . (16) 

I have defined G(r, r') as a free propagator 

V 2 G{r,r') =5(r,r'). (17) 

At boundaries, one sets G(r, r') to be zero. Superscript SI in Eq.16 refers to the first order 
nonadiabatic corrections to the symmetry part of one-particle density matrix. The matrix 
Mf^, or more specifically Mf^(k, r; u, {g v >(T)}, {g v i(T)}) is widely cited in the rest of this 
article. 

Correspondingly, one can also calculate the contribution to the asymmetric component 
of density matrix in the first order adiabatic approximation using Eq.10. Substituting these 
results into the expression for currents, one arrives at 

fdu d 3 k dH a/3 (k,r;{g lf (T)}) d 1 S1 
^ = l^j2^ Dk ^ ®^[^^ Q (k,r;o;,W(r)},W(r)})], 

f duo d 3 k „ dH a i B (k, r: \q n >(T)}) d r 1 , ,<,, , , n , s 

K = \^^f Dk ^- a ^ dr f ;J ^ ^[-^Mg(k,r;^,{^(T)},{^(T)})]. (18) 
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I have neglected a term which does not contribute to the total current because of the 
boundary conditions in Eq.14. In the rest of the article, I will use the notion H a/3 (k,r;T) 
and M^g(k, r; u, T) without showing {g v /(T)}, {g v r(T)} explicitly. 

In the absence of spin-dependent impurity scattering, S z is a good quantum number; 

G^(k, w ) = (u; - e k + a z g^ B B + ir^ 1 )'}- (19) 

The kinetic energy is = h 2 'k 2 /2m — ep and n (ck) is the Fermi distribution of electrons. 

The total charge Q xx and spin (pointing along the direction of Zeeman field or along 
the z-axis) M* x pumped along e x direction per period are evaluated using Eq.18. The final 
results can be expressed in a form similar to Eq.15. 

N N 
»7,£=1,.. »7,f=l," 

n' c ' /s = n'(+) ±n'(-). (20) 

Alternatively one can obtain the charge transport by directly evaluating Eq.15 taking into 
account Eqs. 11,16. 

In Eq.20, I have introduced two antisymmetric tensors, 

/I [To Q Q 

dg v A dg% = - J dT[g v (T)—g ( (T) - g^T)— gri (T)] 

/r rlvdv 1 f) 
J — V ^Tx G{r ^' mr ' )] (21) 

And 

n'(±) = / de k D k ^-[ i- ^^n (e k ± g^ B B)\. (22) 
J de k D k de k 

Here v{e-k) is the one-particle density of states and the volume of structure is V = L x L x L. 
n'(±) depends on the compressibility and n'(±)x^ = n V £ defines the longitudinal adiabatic 
curvatures of spin-up and spin-down electrons. Eqs. 20, 21, 22 are the general results for charge 
and spin pumping in the semiclassical limit. 
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Introducing k F , e F and D as the Fermi momentum, Fermi energy and diffusion constant 
at the Fermi surface respectively, one rescales all quantities in Eq.22, 

|k| = k F a, e k = e(a)e F , D k = h(a)D ; 
g ^ = /(a) g ,n (e k ) = n {e{a)e F ) (23) 

where e(a),h(a) and f(a) are intrinsic functions determined by the energy dispersion, a 
varies from zero to infinity; at fermi surfaces, a = 1 and e(l) = and h(l) = f(l) = 1. 
The longitudinal adiabatic curvatures are 

lc dv(k F ) lc , gn B Bdv{k F ) , 

11 = ^ 77 ' U = a 71 ' 

aei? e F oe F 

^ - - 2 h h < a ^ - /^ISdw (24) 

I am mainly interested in zero temperature results here and in the following sections. If 
e, /, h are assumed to be smooth functions in the vicinity of a = 1 and their dervatives 
at a = 1 are much less than unity, then in the leading order one obtains n lc = 2, and 
tt 1s = d a h(l)/d a e(l). 

Two important general features of spin pumping deserve some emphases. One is that 
the spin pumping current is zero if g v (t) is identical for rj = 1,2, ..N. More generally one 
can show that the pumped charge and spin have to be vanishing if the trajectory of vector 
g = (<7i,<72, •••5'aO in the N-dimensional space encloses a zero area. This is a well known fact 
emphasized in a few previous occasions where charge pumping was studied; the pumping is a 
pure geometric effect determined by two-form curvatures and is absent in a one-dimensional 
parameter space. 

The second feature is that the longitudinal spin current is proportional to the difference 
between adiabatic curvatures (n'(±)) of spin-up and spin-down electrons. Therefore the 
longitudinal spin pumping eflicience is 

Qxx e F ir lc 

In the presence of orbital magnetic fields, one can also evaluate the transverse charge 
and spin pumping current, 
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N N 

(26) 

The transverse adiabatic curvatures which lead to these currents are 

iT c '* s = r n c [u\+) ± n'(-)]. (27) 

Q c = eH ext /mc is defined as the cyclotron frequency of external magnetic fields. Obviously, 
one can also introduce transverse pumping angle 6q and transverse spin pumping angle 9s 
in analogy to the usual Hall angle, 

tan^c = j^- = r f2 c , 

^XXX 



tmes= orr Ton —^ (28) 



Readers can easily confirm these results. In this scheme, the spin pumping current vanishes 
in the absence of Zeeman fields because n'(±) are identical. 

IV. TOPOLOGICAL SPIN PUMPING I: TOPOLOGICAL BEAM SPLITTING 

The key idea of topological spin pumping lies in the fact that the transverse motion of 
electrons is not only affected by usual orbital magnetic fields, or a gradient in Zeeman fields 
but also by spin rotation. Compared to Lorentz forces which act on spin-up and spin-down 
electrons undiscriminately, the topological force induced by spin-rotation does discriminate 
spin-up and spin-down electrons as if they are oppositedly "charged". Indeed as one will 
see, spin-up and spin-down electrons carry opposite charges defined with respect to the 
Pontryagin topological fields. Splitting of spin-up and spin-down electrons in topological 
fields therefore is named as topological Stern-Gerlach splitting (TSGS). So before studying 
the kinetic approach to topological spin pumping, let us offer a qualitative picture of the 
phenomenon of TSGS. 

Apparently, spin rotation doesn't occur in free spaces where the electron spin S 2 is a 
good quantum number. So for TSGS to happen, certain mechanism has to be introduced 
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to rotate spins during transport. It can be achieved by a coupling between electrons and an 
artificial background "magnetic" configuration. To illustrate this idea of TSGS, one studies 
the following Hamiltonian 

H = + Vij) + 9fiBB ° a ' " (r) + Vextir ' T) - (29) 
The unit vector Q, is defined by two angles 6{r) and 0(r) in spherical coordinates; 

Q x (r) = sin6 l (r) cos0(r), 
£l y (r) =sin6 l (r) sin0(r), 

Q z (r) = cos#(r). (30) 

Consider the following coherent states of electrons 

|n(r);+ >= *(r) <g> cos - exp(-i-) exp(tan ^-e^a') | |>, 

2 2 2 

|n(r);- >= *(r)®sin^exp(-i^)exp(-cot^e^(j-)| |> • (31) 

I t> is the spin- up state defined along z axis and o~ is the corresponding lowering operator. 
|f2; ± > are spin-up and spin-down states defined in a local frame where vector e z coincides 
with unit vector fi(r); i.e. 

n(r) -a\ft;± >= ±\Q;± > (32) 

at every point r. Electrons in these states experience X-space spin rotation (XSSR). 

These states are called spin-plus and spin- minus states to be distinguished from spin-up 
and spin-down states defined before. Obviously, plus and minus states discussed above are 
exact eigen states of the local Zeeman coupling and their degeneracy is lifted at a finite B . 

To demonstrate the beam splitting, one evaluates the expectation value of energy oper- 
ator H in these two sets of states. The results in this limit can be conveniently cast into the 
following form 

< U > ± =< tf( r )| J_[ p ± A x (r)] 2 + V(r)|tt(r) > ± gfi B B . (33) 
2m 
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Here p is the momentum operator to be distinguished from the momentum k which is a 
variable. A x -fields are confirmed to be the vector potential of following topological fields 

£\^u = _ £fi\u = —£\vn is an antisymmetric tensor. In both Eqs.33,34,1 use superscript X 
to refer to X-space gauge fields. And in Eq.33, I have neglected terms which are identical 
for |± > states. This form of the kinetic energy of spin-rotating electrons was previously 
derived to demonstrate interactions between quasi-particles and spin fluctuations in triplet 
superconductors [18]. 

Therefore the total forces exserted on spin-plus and spin- minus electrons are 

F±= <*j-^W ±T *xP±^W|*>. (35) 
or in 

The last identity holds when the spin rotation is adiabatic so that transitions between 
Zeeman split spin-plus and spin- minus states are negligible. I.E. 

; ' °' >(r ) . / m **(r)|-*(r)| « g^ B B . (36) 



2m dv dv 

One can verify that when this adiabaticity condition is satisfied, a • n is an approximate 
good quantum number. 

In the semiclassical approximation employed below, I further assume that 

|/m*»^(r)!»|^|,|^| (37, 

So over the wave length of electrons, spin rotation is negligible, 

In addition to a term proportional to the field gradient of scalar fields, there is a new 
force perpendicular to the velocity of electrons similar to the Lorentz force. More important 
as indicated in Eq.33, spin-plus and spin- minus electrons carry opposite topological charges; 
therefore, the corresponding forces are in fact along exactly opposite directions. This shows 
that an XSSR does affect the orbital motion of electrons and does differentiate spin-plus 
states from spin-minus states. It therefore leads to the promised phenomenon of TSGS. 
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It is important to further emphasize here that to observe TSGS, the background configu- 
ration has to be topological nontrivial (see more in the next section) so that T x is nonzero. 
To highlight the relevance of topology of spin configurations, let us consider spin states 
defined in Eq.31 where fi(r) corresponds to a hedgehog configuration, 

0(r) = j^. (38) 

Our calculations show that forces acting on two spin-rotating electrons given in Eq.31 are 
equivalent to forces exserted on two oppositely charged particles in a resultant magnetic 
monopole field 

(39) 



2|r| 3 

Before leaving this section, I generalize the argument to momentum space topological 
effects. I then consider two orthogonal k-space wave packets; spins in these states are 
pointing at either k or — k direction. Let us define 6(k) and 0(k) as 

n x (k) = sin#(k) cos0(k), 
n„(k) = sin#(k) sin0(k), 

fi*(k) = cos#(k). (40) 

and fl (k) = k/|k| is a unit vector along k. 

|n(k);+ >= ^(k)cos-exp(-i-)exp(tan^e^)| |>, 

2 2 2 

|fi(k);->=*(k)exp(-Asiii^exp(-cot^e < *0| T> • (41) 

As before I assume these spin-plus and spin-minus states are split by an effective Zeeman 
field B . 

An electron in such a wave packet experiences k-space spin rotation (KSSR), or spin 
rotation that depends on its momentum. For the same reason mentioned before, I further 
assume adiabaticity in spin rotation and use a semiclassical approximation. This requires 
that 
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gH B B ^>\—(j>(k)-—H(r)\, 

it r*n \ d t/1 M ,^( k ), ,^( k ) 

|/^(k)-^(k)|»|^|,|^ ; 



(42) 



The group velocities of these wave packets are 



v± = % < [r, W(r,p)]_ > ± =< *(k)|^ =F [V x A*(k)] x ^ ;r) ^(k) > . (43) 

Superscript K is introduced to specify the k-space gauge fields. So the velocity does acquire 
an additional nontrivial transverse term in the presence of an external field gradient and 
topologically nontrivial fields T^(k) of which A K {k) is the vector potential. A calculation 
shows that 

- 4^«(k) • x _. (44) 

And more important the spin-up and spin-down electrons defined in the local frames drift 
in an opposite direction once again leading to TSGS. 

The general property of electrons illustrated in Eq.43 has been noticed in a few different 
occasions. In Ref. [9,10], T K fields were expressed as Chern- number density of electron 
states which was first introduced for the studies of quantum Hall and fractional quantum 
quantum Hall conductances [6]. In Ref. [14], T K fields are from a topological monopole in 
the k-space. 

The general form of topological fields obtained in Eq.44 is a generalization of standard 
Pontryagin fields defined in the X-space. One easily recognizes that topological fields dis- 
cussed here are equivalent to Berry's two-form fields defined in an external parameter space 
[4]. They generally represent holonomy of parallelly transporting eigen vectors in the Hilbert 
space [5] , in this particular case the holonomy of transporting a spin-up or spin-down eigen- 
vector defined in local frames. I will discuss topological spin pumping in the presence of this 
TSGS in section VI. 



15 



V. TOPOLOGICAL SPIN PUMPING II: ADIABATIC SPIN TRANSPORT IN 

THE PRESENCE OF XSSR 



I now take into the kinetics which leads to spin rotation while electron wave packets 
propagate in the X-space. I study the spin pumping phenomena in this case via applying 
kinetic equations derived in section II. 

Consider electrons coupled to a background spin configuration or an artificial magnetic 
field with uniform magnitude but with spatially varying orientation. One models electrons 
with the following Hamiltonian 

if (k, t ,T) = ^-hf- gpL B [B n(r) ■ a + Ba z \ + V ext (r, T). (45) 

again f2(r) is a unit vector representing the orientation of an internal exchange field in metals. 
The scattering of impurity potentials is taken into account in elastic collision integrals in 
Eq.4. I assume g\iBB§,g\i B B are much smaller than the fermi energy ep. 

To facilitate discussions, one introduces the following unit vector n(r) which defines the 
direction of net Zeeman fields in the above equation 



+ og + (n, + if 



n 



fL + f 

n 2 = . (46) 

where I = B/B . 

Alternatively, similar to Eq.30, in spherical coordinates one introduces the following 
characterization of n, 

n x = sin^(r) cos0(r), 
= sin^r) sin0(r), 
n z = cos^(r). (47) 
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I use superscript tilde to distinguish the spherical coordinates 9,4> for n from 6>,0 for fi. 
Eq.46 indicates that 

t i / \ n sin#(r) , . 

= r ),6 = arctan \> 48 

cos#(r) + I 

Introduce a local spin rotation such that 

U- 1 {r){n{v)-a)U{r) = (j z . (49) 

The Hamiltonian becomes 

H(k,v,T) = (k ~fJ (r))2 -VF- 9^bB [Q 2 x + Q 2 y + (SI, + J) 2 ] 1/2 ^ + K,*(r,T). (50) 

Here SU (2) gauge fields generated by pure spin rotation are 

A?=^- 1 (r)J-U(r)=Af^. (51) 

7 = x,y, z. To simplify the formula, in this equation and the rest of article I do not show 
spin indices explicitly. A direct calculation yields 



A x y 



ld9(r) 



" 2 9r M ' 



Af = -oos6l(r)-^; (52) 



And full covariant SU(2) fields vanish as one should expect for pure spin rotation. 



To proceed further, one notes that the degeneracy between spin-up and spin-down states 
in the rotated basis or spin-plus and spin- minus states is completely lifted by various Zeeman 
fields. One again assumes that spin rotation is slow in the X-space so that the adiabaticity 
specified in Eq.36 can be satisfied. In the adiabatic approximation, one neglects transitions 
between spin-plus and spin- minus states and set off-diagonal gauge potentials A Xx ,A Xy as 
zero. Therefore in Eq.50, one only keeps A Xz which yields the usual Berry curvatures for 
spin-plus and spin- minus states. Corresponding U(l) gauge fields are 
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So in this limit, only the ^-component of SU(2) gauge potentials survives to contribute 
to pumping currents; it also defines well known Pontryagin type U(l) fields 



jx-^v S ** 1 cind(r)\ d§ir)d ^ r) ^>W(r) f 54 ) 



6* and 4> again are two spherical angles of n(r) in Eq.46. 

In the rotated basis, the structure of equation for one-particle density matrix should be 
identical to the one in the previous section. However according to a general consideration in 
the semiclassical transport theory, the following transformation takes place in the equation 
of motion for electrons in the presence of XSSR. 

dH(k, r; T) 



dr dt 

^^^M(M;t)|., (55, 

is the usual Poisson bracket defined with respect to canonical coordinates {u,v}. k is 
the electron momentum 

k = p + A Xz (r)a z ; (56) 

and p and r are a pair of canonical coordinates. In the semiclassical approximation, the 
new equation for the one-particle density matrix therefore looks as 

dp(k,r;u;,T) n dH(k,v;T) x dH(k,r;T)^ d 



dH(k,r;T) d dH(k,r;T) d , „ _ ~ T „ 

V ' ' ; _ + \' ' ® P(k, r; T) = l CI p(k, r; T) (57) 



<9k M <9r M <9T &u 



Following Appendix A, the transverse pumping currents are given as 

J » -I I ^ mr{ ^ 8 ^7^^ [ ^ MS1 (k ' * r)1} 
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where tr is only taken over the spin space. This is the central result for topological transverse 
spin and charge pumping in the presence of XSSR. 

Let me now again consider external perturbations specified in Eq.13. To address spin 
pumping, I consider a background spin configuration of square half-skyrmion lattice given 
below, 

/ \ • / \ TT Z — ZUl,^) . /™\ 

n x (r) + tn y (r) = [[ -== =;z = r x + tr y . (59) 

h,h yj\z - z(l 1: l2)\ 2 + y s 



z is introduced as a coordinate in the 2D plane. z(l±, I2) = ha + il 2 a represents a lattice site 
with Zi )2 as integers and a the lattice constant. I also assume that As <C a. The average 
topological fields of this lattice are 



7T 



Eif = £0 = (60) 



The adiabaticity condition in Eq.36 requires that 



< gfi B min{B ,B}. (61) 



The equilibrium density matrix is still given in Eq.19. In the longitudinal direction the 
results are the same as in Eq.24 and I will not repeat here. Furthermore, some straightfor- 
ward calculations lead to the following expression for total charge and spin transported in 
a transverse direction y per period 

N N 

Q yx = n tc £ xA,M* x = n ts £ xA; 

»7,£=1,.. r),£=l,.. 

n tc ' ts = r n ct [ii l (+) T n'(-)], n ct = (62) 

m 

And Qtf is introduced as an effective cyclotron frequency of topological fields S . 

The transverse charge pumping angle 9c and transverse spin pumping angle 9s are 

tan 9 C = r il ct F{-5-), 

6f d 

t&n9s = T Q ct . (63) 
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F(P) is a function of (3\ it approaches unity as (3 becomes much less than one. The transversal 
spin pumping efficiency in this case is 

e t = = tangg = e F F -i,-go 
Q yx tan 8 C g^sB B 

It is important to notice that e* diverges as the external Zeeman field B goes to zero 
signifying zero transverse charge pumping. In practical cases, the topological spin pumping 
is always accompanied by small transverse charge pumping current, a unique feature of 
topological fields. This also occurs in the scheme discussed in the next section. It is however 
in contrast to spin pumping of polarized electrons discussed in section II where the transverse 
spin pumping current is negligible compared to the transversal charge pumping current. 
This feature originates from the fact that spin-plus and spin- minus electrons carry opposite 
topological charges. In the absence of polarization, there are equal amplitude of currents of 
spin-plus and spin minus electron flow in opposite transverse directions as a result of TSGS. 
So the total charge pumping is zero but spin pumping current flows. 

In all topological pumps discussed here and below, I have found that spins are pumped 
out by applied a. c.gate voltages because of a spin configuration which yields either nonzero 
as in Eq.58 or nonzero E^ as in Eq.85. I intend to call these topological configurations 
as topological motors in spin pumps. 

VI. TOPOLOGICAL SPIN PUMPING III: ADIABATIC SPIN TRANSPORT IN 

THE PRESENCE OF KSSR 

To illustrate the effect of KSSR, I start with a 3D toy Model and discuss the topological 
mechanism. In the second half of this section, I study the topological mechanism in more 
realistic models for electrons in semiconductors. 

In all subsections here I assume that spin-orbit splitting and Zeeman field splitting are 
much smaller than the fermi energy but can be comparable between themselves. Zeeman 
splitting is introduced to ensure that spin degeneracies at k = are lifted and the adiabatic- 



■ 



(64) 
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ity holds for every state below fermi surfaces. The zero field limit should be taken when the 
adiabaticity conditions in Eq.99 are satisfied. 

In the presence of impurity scattering, to ensure adiabaticity I always assume that the 
impurity potentials are weak compared with the splitting between two spin bands either due 
to spin-orbit coupling or due to Zeeman field splitting. I neglect therefore the interband 
transitions due to nonadiabatic corrections. Because of this reason, I only consider the limit 
of strong spin-orbit coupling and present results at zero temperature. 



In this section I consider electron spins coupled to the momenta of electrons and spin 
rotation occurs when a wave packet propagates in the momentum space. A propagating 
wave packet in the k-space corresponds to an accelerated electron. The artificial model 
I introduce here to study topological spin pumping can be considered as a mathematical 
generalization of the Luttinger Hamiltonian [19] to spin- 1/2 electrons. 

Consider the Hamiltonian 



In the toy model there are two spin dependent terms; the term proportional to <r z represents 
the Zeeman splitting of electrons in the presence of fields along the z-axis and f2 • a term 
characterizes a collinear spin-orbit correlation. 

To understand the topology of spin configurations in Fermi seas, I again introduce a unit 
vector n(k) as a function of f2(k) defined in Eq.46 in the previous section. Especially, 



A. A toy Model 



if (k, r; T) = e k - ^ BJ B (|k|)[/(|k|)a 2 + fi(k) • a] + V^(r; T); 



(65) 



I and Bq are functions of |k| and f2(k) is a unit vector along the direction of k; 




(66) 



0(k) = 0(k), 



6 = arctan 



sin#(k) 



(67) 
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Again 6>(k), 0(k) are two spherical coordinates of unit vect il(k) = k/|k| and differ from two 
angles of n(k), 9(k) and 0(k) (see Eq.30). 

One then considers a configuration of n(k) on a sphere at a very large momentum which 
naturally defines a mapping from an external S 2 sphere in the k-space to a target S 2 space 
where n(k) lives. The topology of electron spin states in Fermi seas is therefore characterized 
by ^(S* 2 ), the second group of target space S 2 . 

Let us further introduce T K fields defined in Eq.44 with f2(k) replaced with n(k). The 
winding number of a mapping or configuration can be characterized by the flux of T K fields 
through a large surface. It is easy to verify that 



where the surface integral is taken over at |k| — > +00. This shows that spins of electrons 
form a monopole structure. This is not surprising because at very large k or small I, n(k) 
is identical to fi(k) which always points outward along the radius direction. 
Again introduce a k-space local spin rotation such that 



H(k, r; T) = e k - gfi B B [Q 2 x + Q 2 y + (Q z + I( Pk )f] 1/2 a z + V ext (r - A K ; T). (70) 




(68) 



f/- 1 (k)n(k) -aU(k) 



(69) 



Under the spin rotation, the Hamiltonian becomes 



As before, SU(2) gauge fields are generated under the spin rotation 




(k)= l U- 1 (k)—U(k) = A^a<; 



(71) 



in a fixed gauge, one obtains 





1 d9(k) 



(72) 
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And the full SU(2) fields again vanish. 

In linear responses, the effective Zeeman field splitting between spin-plus and spin- minus 
states in Eq.70 is stronger than external perturbations. Further, I require that the energy 
splitting is also stronger than impurity potentials. So the adiabaticity condition in Eq.42 is 
always satisfied. I therefore set A Kx , A ky to be zero again to neglect transitions between spin- 
plus and spin- minus states. In Eq.70 I only keep the ^-component of SU(2) potentials A^f 
which yields to Berry's phases of plus and minus states. The corresponding z-component of 
reduced SU(2) gauge potentials which enters our results below is 



(k) = — ^ " (73) 



I note that in the presence of spin-orbit coupling, impurity scattering combined with 
A Kx,Ky components of SU(2) gauge fields does lead to transitions between different spin 
bands. The adiabaticity condition in this case however is sufficient to ensure that these 
contributions are negligible. So the Born-approximation employed in this article is valid in 
the strong spin-orbit coupling and finite Zeeman field limit where the splitting between plus 
and minus spin bands at any momentum k, 

2^ B r|k|[^ + n 2 y + (n, + /(p k )) 2 ] 1/2 (74) 

is much stronger than the impurity potentials. 

Where the adiabaticity condition is satisfied, one shows that 



^ v ^ 2 v ' 9k M dk u 



(75) 



In spherical coordinates where k = (p k , 9, 0), one has the following explicit results 

1 + / cos 9 1 



l + / 2 + 2cos^/) 3 / 2 2p 2 k 

yKz = 1 d/(p k ) 

p4> - (1 + /2 + 2 cos 91 fl 2 2p k dp k 



and I is a function of p k (= |k|) as defined in Eq.66. 
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It is convenient to redefine topological fields in terms of T^f, 

T?(k) = ^E**(k). (77) 

One finds the following asymptotic behaviors at large and small momenta, 

K[l-^cos0]e p - Hsm9e e , when p k > A m ; 



T K (k) 
And here 



r\ V -L Kj\J>3 W n 3 

2p£ L Pk J p 2p£ 

iyk-e z + ^-§-[(1 — 3cos 2 6')e p + 3sin#cos#e e ], when p k <C A r 

ZA™ ZA™ 1 



(78) 



A m = ^ (79) 



defines the core of anisotropic monopoles when the Zeeman fields are present. 

When the Zeeman field B is set to zero or I = 0, unit vector n(k) coincides with Q(k) 
and 4> = 0(k), 9 = 0(k); therefore, S^f vanishes. Eq.76 in this limit indicates familiar 
isotropic monopole fields in the k-space. In the presence of finite Zeeman fields, topological 
fields are not strictly isotropic because the inversion symmetry z — > — z is broken by external 
Zeeman fields. Topological fields are along the z-axis at small momentum limit but approach 
isotropic monopole fields at large momenta. The crossover takes place at \ m . 

Under KSSR, the following transformation should occur in the equation for the density 
matrix, 

8H(k, r; T) <9k 



dr dt 
dH(k, r; T) dr 



dk dt 
The electron coordinate in the presence of KSSR is 



= [k, J ff(k,r;T)] X)k ; 

[r,tf(k,r;T)] Xik . (80) 



r = x - A Kz (k)a z ; (81) 



x, k are a pair of canonical coordinates in this case. 
The corresponding kinetic equation becomes 

9 ^^ + ^) + ^.w- BE ^ r ^ B 



dT ^ v ' dr v J <9r M 



dH(k,r;T) d dH(k,r;T) 9. ,, ^ „ T 

+ ^^^}®p(k,r;cu,T)=J a/ p(k,r;a;,T) (82) 



dr^ dk^ dT du~ 
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The charge and spin current expressions transform accordingly; in the rotated basis one 

has 

J(r ' T) = J t (2^ Tr{K ^ (k) + S - (k) ® dHi dC T) ] P(k ' r; T)} 

JZ(r ' T) = / ^ Tr{n(k) ' ' kn " (k) + * (k) ® ^7^ ] ^ r; r)} (83) 

Eq.83 can be used to analyze contributions to the spin and charge pumping currents from 
different part of Fermi surfaces. 

Let us define plus and minus Fermi seas as shown schematically in Fig. 2. In the plus 
Fermi sea, electron spins are along the direction of unit vector n(k) and in the minus Fermi 
sea spins are along the opposite direction of n(k). The corresponding Fermi surfaces are 
named as the plus and minus fermi surfaces. For a system with an inversion symmetry each 
fermi sea has zero overall polarization when the Zeeman field is absent. In the rotated basis, 
these two fermi seas correspond to spin-up and spin-down ones. 

Consider electrons subject to pumping potential gradient along the negative x-axis. One 
easily finds that spin-up electrons at the north pole of plus Fermi surface are subject to a 
drift along the y-direction while spin-down electrons at the south pole of plus fermi surface 
are subject to a drift along the minus y-direction. For the same reason electrons at the 
north and south pole of minus Fermi surface are subject to a drift along minus and plus y 
direction respectively (see Fig. 4c). 

Following these results one also finds that in the absence of spin polarization, charge 
pumping currents carried by spin-up electrons at either the north pole of plus Fermi surface 
or the south pole of minus Fermi surface flow in the exactly opposite direction of charge 
pumping currents carried by spin-down electrons in the south pole of plus Fermi surface 
and the north pole of minus Fermi surface. So while the spin pumping current flows, the 
net charge pumping vanishes. Only when electrons are polarized or B is nonzero, charge 
pumping is possible along the transversal direction. 

Given the expressions for p A in Eq.ll and the current expressions in appendix B, one 
evaluates the transverse charge and spin pumping currents, 
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W ^/s5^»' k » ' ^« * d -^^M^,T)Y (85) 

Taking into account Eq.75, one then arrives at expressions for charge and spin pumping 
currents. The longitudinal spin and charge pumping are still given by Eq.20. The transverse 
spin and charge pumping currents are more involved. To evaluate the spin and charge 
current, one notices the following identities 

Tr(^) = 2E**, TrKE*) = Tr(^) = (86) 

according to Eq.75. Final expressions for spin and charge pumping currents only depend on 
the z- component of reduced SU(2) fields. For this reason, one is able to obtain a rather 
simple form for transversal adiabatic curvatures defined in Eq.26 

= W + * + («■ + / (87) 

The topological motor S^f in this case is an anisotropic monopole discussed in Eq.75. I 
only present results in the following weak Zeeman field and strong Zeeman field limits. 

Since the topological fields have distinct large momentum and small momentum asymp- 
totics, the topological pumping has strong dependence on Zeeman fields. The spin and 
charge topologically pumped out per period are again given by Eq.26; the transverse adia- 
batic curvatures H tc ' ts are calculated and results are 

n*. = 7r f(**)i_J_iifc ; 

K \Jn lc D m 



U lc , n lc are provided in section III; D is the diffusion constant and m is the electron mass. 
I should mention that when the Zeeman field vanishes, IT' C goes to zero but Tl ts remains 
finite; as pointed out before, this is a distinct feature of a topological spin pump where an 
electron beam splits because of TSGS. 
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B <C Tk F This corresponds to a limit where the core size of monopole \ m is much smaller 
than k F , 

^) = ^ln^,vr-(^) = o(^). (89) 

A m 1Z A m A m K F 

B 3> Tk F This corresponds to a limit where the core size of monopole \ m is much larger 
than hp, 

A m Z A m *i /\ m 

In deriving these results for 7r' s ' ic , I have neglected the k-dependence in diffusion constants 
and derivatives of density of states. 

The corresponding charge and spin pumping angles are 



tan#c* = n 



tc ,k F 1 gfi B B 1 



tan0 5 = 7r' s (^)^--^— . (91) 



\ m 7T lc e F D^va 
k F ^l 1 
\ m tt 1c D m 

Finally, the transversal spin pumping efficiency is 

r ts( kp 



, tan9 s 7T e F 

(92) 



tanOc vr* c (^) gf jL B B' 



B. The 2D Rashba Model 

In the previous section, I discuss the topological spin pumping due to the k-space spin 
rotation in an artificial model. Now I turn to more realistic models for semiconductors. And 
I limit ourselves to 2D cases. Spin-orbit coupling can be either due to the Dresselhaus term or 
the Rashba term [20-23]. In the later case, or in the Rashba model for 2D semiconductors, 
the spin-orbit coupling has a particularly simple form because of either a bulk inversion 
asymmetry or a structure inversion asymmetry [24]. I start with discussions about this 
model. 

In the 2D Rashba model, the spin dependent Hamiltonian can be presented as 
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H R = - gf x B B (k)[I(k)a z + e z ■ Q x a}. (93) 

In bracket [], the first term is due to Zeeman fields and the second one is the Rashba coupling 
term. As in Eq.65, I have defined J(k) = B/B (\\c\) and B — T|k|. f2(k) is a 2D unit vector 
along the direction of k, 

k 

Vl x = cos0(k) = j^j, 

^ = sin0(k) = -^. (94) 

Here 0(k) is a polar angle of k vector in the 2D plane. 

To characterize the spin configurations in (k x , k y ) plane, I study the unit vector n(k) 
defined as 

n = — ■ 

x VT+P' 



n y 



n z = - 1 (95) 
One also obtains simple results for spherical angles 0(k),<fi(k) of n(k) in this case, 

= 0(k) + ^, 

= arctan — . (96) 
^ (1^1) 

One notices that at k = 0, I(k) becomes infinity; as a result unit vector n(k) points at 
the direction of e z because of Zeeman fields. At the large k limit 6 approaches tc/2; n relaxes 
and lies in the equator plane of two sphere S 2 . So 9 varies from to 7r/2 as one moves away 
from the center of Fermi seas while = 0. This behavior of unit vector n implies a meron or 
half-skyrmion in the 2D momentum space. Merons have been proved to play important roles 
in Yang-Mills theory as well as in quantum magnetism [25,26]. The size of half-skyrmion 
outside which the spin polarization along ^-direction becomes unsubstantial is 

A, = f • (97) 
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To confirm the peculiar topology of electron spin states one examines the homotopy class 
of mapping from (k x , k y ) space to a two sphere S 2 defined by n(k). Consider T K defined in 
Eq.44 in terms of n(k) vector instead of f2. The winding number can be easily calculated 

W = ^- [dSe z -T K = l r /2 d0sin0 = i (98) 
Z7T J 2 Jo 2 

which precisely shows a meron in k x — k y plane. 

In general, T is a quantity which can be controlled by an electric field [23]. X s is a function 

of both Zeeman fields and applied external gate voltages which offers great opportunities 

to manipulate the meron structure and control topological spin pumps discussed below. In 

this model, the adiabaticity condition in Eq.42 requires that for each k, 



dV ext (r) 1 1 ,dV im (r) 1 1 

l?™n{ — , TTT/' I — a \min{ — , — } -C gfji B max{B (\k\), B). 



dr \ s |k| dr A s |k| 



(99) 



Here V im is the impurity optential. The sufficient condition for Eq.99 to hold is that 



M ext dV im (gy B B) 2 
dr dr T ' 1 } 

Furthermore, the size of systems has to be larger than min{A7 1 , |k| -1 } to ensure the semi- 
classical approximation [28]. 

At last I would like to emphasize one more time that the impurity potential has to be 
weak compared with either the Zeeman splitting or the splitting due to spin-orbit coupling so 
that the adiabaticity condition can be satisfied. Therefore, the transitions between the split 
spin bands are also negligible in the adiabatic limit. As mentioned in a few occasions in this 
article, for this reason the Born approximation is always valid in the adiabatic approximation 
[29]. 

Similar to the procedure introduced in the previous section, it is possible to introduce spin 
rotation to diagonalize this Hamiltonian. In the rotated basis, spin-up and spin-down states 
are split by a combined field of external Zeeman splitting and internal spin-orbit coupling 



with the following strength, gusBoy 1 + J 2 (|k|). The resultant U(l) gauge fields in the 
two-dimension k-space are vortex-like. I present the result in polar coordinates k = (p k , 0); 
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A?2(k) " 2p k VT+P : 



The asymptotics for the vector field T K are 

-e 2 , when p k > \ s ; 



T (k) = { 



k 



(102) 

_ , ^ Id wVion n. <^ \ 

2A?^ 2A? 



1 (1 - §r£)e z , when p k < A,. 



In the absence of Zeeman fields, I find 9 = 9 and S Xz is zero everywhere in the momentum 
space except at k = 0. That is 

K z = #)■ ( 103 ) 

Note that the topological fields though zero everywhere are singular at the origin of the k- 
space. In general, topological fields are negligible when p k is much larger than the skyrmion 
size A s . However, I find this is sufficient to produce a spin pumping current even if the 
skyrmion size \ s is much smaller than the Fermi radius hp- The topological motor in this 
model is a positive meron (see Fig. 3). 

Let us again define plus and minus Fermi seas. In the plus Fermi sea, all electron spins 
are along n(k) while in the minus Fermi sea all electron spins are along — n(k) direction. 
In the rotated basis, the plus and minus fermi seas become spin-up and spin-down fermi 
seas respectively. In plus and minus Fermi seas, electrons subject to external pumping fields 
along x-axis drift along plus and minus y-direction respectively, as shown in Fig.4 (d). 

The spin and charge topologically pumped out per period are again given by Eq.26; the 
transverse adiabatic curvatures H tc ' ts are given in the following equations 



2vr D k 

I have used subscript R to refer to the adiabatic curvatures in the Rashba model. Taking 
into account the profile of topological fields in Eq.101, I obtain the following results for the 
transverse adiabatic curvatures 
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n£ = 7 r*(i*)i--J_n ,c ; 



n£ = 7r te (^)^-^^V. (105) 



(see section III for discussions on Tl lc ,ir lc ). ir ts and 7r te are calculated using rescaled param- 
eters introduced in Eq.23, 

77 {b) = 2J da ^TW)a{l + Pf Oa n °' 

^( 6)= /daMl 1 J( ^y ^. (106) 
v 7 y /i(a) al + P da da v 7 

I again present results in the limit of strong and weak Zeeman fields. 
B <C Tk F This corresponds to a limit where the core size of meron A s is much smaller 
than kp. 

B 3> TkF This corresponds to a limit where the core size of meron X s is much larger 
than &p. 

vr-(^) = = ^) 2 - (108) 

The corresponding charge and spin pumping angles, the transversal spin pumping effi- 
ciency are still given by Eq.91, 92; for the Rashba model, ir ts and ir tc calculated above should 
be used to determine the angles and efficiency. Without losing generality, I again have ne- 
glected the k-dependence in Dk and <9z/(ek) /de-^ in deriving results for 7r ts ' tc in this section; in 
the 2D model, I further choose to work in a limit where dv(ey s ) /de^ is nonvanishing because 
of band structures so that the longitudinal charge pumping is nonzero. 



C. The 2D Dresselhaus Model 



In the Dresselhaus model, the spin-orbit coupling and Zeeman coupling are given as 

H D = - gf i B B [a~n(k) + a z I(\k\)}; (109) 
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here A~B is defined as A X B X — A y B y . Again one introduces a unit vector n(k) to specify 
spin configurations in fermi seas. The unit vector n(k) characterized by 9(k) and </>(k) is 
given in terms of /(|k|) and in the following equations, 

= -0(k), 

^ = arctan 7(M' (U0) 

Calculations for spin pumping current are identical to those in the previous section. The 
winding number of the configuration defined by n(k) in this case is 

w=i LJ dSez ' TK = ~\ \T d§s[n§ = ~\ (ni) 

representing a meron with negative one half skyrmion charge. This is topologically distinct 
from the spin configuration in the Rashba model. Naturally, all topological fields in this 
model are pointing in the negative z-axis. 

In the Rashba model I find that the topological motor is a positive meron in fermi seas 
while in the Dresselhaus model the motor is a negative meron. I anticipate that spins should 
be pumped out in an opposite transverse direction in these two models. This is as well true 
for charge pumping currents in two limits when a Zeeman field is present. So both topological 
spin and charge pumping currents flow in an opposite transverse direction compared with 
currents in the Rashba model. 

More specifically, I find that transverse adiabatic curvatures in the Dresselhaus Model 
(I use subscript D to refer to this model) are related to those in the Rashba model (I use 
subscript R for this model) via the following identity 

u ts,tc = _ U t£tc ( U2 ) 

This is an exact result as far as the adiabaticity conditions in Eq.36,42 are satisfied and is 
independent of the spin-orbit coupling strength in the 2D Rashba and Dresselhaus models. 

In all models, I have found that topological spin and charge pumping are suppressed by 
strong Zeeman fields because of spin polarization. This is a general feature of topological 
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pumps; the topological fields are absent when electrons are completely Zeeman polarized. 
Furthermore, the topological charge pumping has a maximum when the Zeeman field is 
comparable to spin-orbit fields, i.e. B ~ Tkp. In the appendix C, I have found similar 
effects on the spin Hall and Hall conductivity. In the case where both the Rashba term 
and Dresselhaus term are present and controllable, it is interesting to understand how a 
topological pump reverses the direction of its spin flow. 

VII. COLLECTIVE RESPONSES OF VARIOUS FERMI SEAS TO PUMPING 
POTENTIALS: TRANSVERSAL DISPLACEMENT VERSUS SELF-TWIST 

I also would like to emphasize that different schemes of spin pumps discussed in this 
article correspond to different responses of Fermi seas to adiabatic perturbations. To study 
displacement of Fermi seas or more general deformations of various Fermi seas, one examines 
the adiabatic displacement when pumping fields are applied along the x-axis. It is useful 
to introduce the displacement vector-tensor d(k) to characterize collective motion of fermi 
seas, 



In the adaibatic approximation employed in this article, only the diagional part of dis- 
placement tensor is nonzero. One studies these elements to analyze the responses of Fermi 
seas to pumping fields. For weakly polarized electrons in orbital magnetic fields, spin-up and 
spin-down Fermi surfaces are displaced in the same longitudinal and transverse directions; 
electrons in each Fermi sea drift with almost the same velocity when the Zeeman splitting 
is much smaller than the Fermi energy. Following the calculations in section II, 



5p(k,r;T) =mvd(k)^-^, 
5p(k,r;T)= |^[p(k,r;c,T) 




p (k;c)],p (k)=/^p°(k;c) 



(113) 




(114) 



The responses of Fermi surfaces are summarized in FIG. 4. 
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In an XSSR based pump, the plus and minus fermi surfaces are displaced along the same 
longitudinal direction but along opposite transverse ^-directions. 

, To dV ext (r, T) 2 r o v x 21 , 11C .\ 

d " = "m dr v 5 ^ + a ( 115 ) 

For KSSR based pumping, Fermi seas respond in very distinct ways. In the generalized 
Luttinger model, the Rashba model and the Dresselhaus model, Fermi seas only experience 
displacement in the longitudinal direction; 

T Q dV ext {r,T) 

^ = ~m dr v 5 ^ x (U6) 

And Fermi seas are not displaced along the transverse t^-direction. In this regard, there is 
a fundamental difference between spin pumping of polarized electrons, an XSSR spin pump 
and a KSSR spin pump. In the later case, the one-particle density matrix, surprisingly, 
doesn't develop an asymmetric component along the transverse direction. The spin current 
therefore has characteristic of persistent currents which doesn't involve distortion of Fermi 
seas. This observation appears to be consistent with a proposed analogy between a spin 
injection current in the Luttinger Model and a supercurrent in Ref. [15]. 

However, in this case the group velocity of electrons acquires a transverse term in the 
rotated basis along the transverse ^-direction; the dispersion of group velocity tensor is 

v, = t,n fcM + (7*E£(k)*^. (117) 

Consider the toy model in section VI A, in terms of group velocities one finds that the upper 
half of the plus Fermi sea is displaced along the positive v ^-direction and the lower half twists 
along the minus tydirection. For the minus Fermi sea the upper half twists along minus v y 
direction and the lower half twists along the plus i^-direction. Therefore plus, minus Fermi 
surfaces experience two distinct self-twists. All these occur while there is no displacement 
of fermi seas along the i^-direction. 

This is as well true for the Rashba model and Dresselhaus model. The displacement 
vector tensor doesn't have a transverse component; only the group velocities of electrons 
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develop a t^-component. Following Eq.117, the dispersion of group velocity is given by the 
calculated meron fields T K in section VI C. 

So I find that in an XSSR spin pump, Fermi seas experience periodical displacement 
along the ^-direction and at any moment, spin-plus and spin- minus fermi seas are displaced 
towards two opposite points in the v y axis. In the KSSR pumping scheme on the other 
hand, Fermi seas do not have periodical transverse displacement; however group velocities 
are renormalized by external fields. These responses of Fermi seas lead to topological spin 
pumping phenomena. 

VIII. CONCLUSION 

To summarize, I have developed a kinetic approach to analyze topological spin pumps 
in some details. In all topological spin pumps discussed here, spins as well as charges are 
pumped out in a transversal direction by various topological spin configurations which I have 
called topological motors. I have introduced the spin pumping efficiency and spin pumping 
angle to characterize topological spin pumps which could operate in the absence of spin 
polarization through a mechanism of TSGS. 

Following Eqs.63,91 for a given longitudinal charge pumping current, the transversal spin 
pumping angle for XSSR based pumps 9s is bigger in clean systems than in dirty systems; 
this is similar to the usual Hall angle in dirty metals. Transversal spin pumping in this case 
is attributed to distinct displacement of spin-plus and spin- down fermi seas. On the other 
hand, following Eq.91 for KSSR based pumping, the transversal spin pumping angle 9s is 
bigger in dirty structures. Transversal topological spin pumping in this case occurs when 
Fermi seas are not displaced along the transverse direction. 

A unique feature of topological pumps is its Zeeman field dependence. We have illustrated 
that when Zeeman fields are much stronger than spin-orbit coupling fields, the topological 
pumping mechanism is strongly suppressed. In the 2D Rashba model, a topological spin 
pump is driven by a meron living in Fermi seas while in the Dresselhaus model, spins are 
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pumped out in a transversal direction by a meron with a negative half skyrmion charge. The 
suppression occurs as the size of meron increases. This general feature also appears in the 
spin Hall and Hall currents in these models. In a subsequent paper, I am going to discuss a 
design which can be used to establish a topological spin pump in laboratories. 

Finally I would like to point out that generalization to various mesoscopic quantum 
systems is possible. A mesoscopic mechanism of charge pumping was proposed and studied 
in a few recent works [27,30,17]. The possibility of pumping charges in an open mesoscopic 
structure was first noticed in Ref. [27]. In an extreme quantum limit, external potentials 
can be applied to manipulate coherent wave packets rather than to vary chemical potentials; 
the resultant pumping current therefore is greatly enhanced at low temperature limit and 
decays rapidly at a few hundred mk to one kelvin range. Such a quantum electron pump 
has been achieved in a remarkable experiment [32]. 

A few interesting aspects of coherent charge pumping were addressed in later works. The 
symmetry of charge pumping in the quantum limit was studied in Ref. [33]. The issue of 
counting statistics was raised and addressed in a fascinating work [34] ; the counting statistics 
might shed light on the issue of dissipation in an adiabatic pump and remains to be studied 
in experiments. Pumping of coherent Andreev states was analyzed in Ref. [31]. At last I 
should also mention that various other proposals for spin pumps have been made by different 
groups; see for example [35,36,38,37,39,40]. I would like to refer readers to those works for 
detailed discussions in various limits. In chaotic quantum dots where current experiments 
are carried out, an adiabatic spin pump of coherent polarized electrons was proposed in Ref. 
[36]; the effect of spin-orbit scattering has been further studied in a recent article [40]. 
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APPENDIX A: CURRENT EXPRESSIONS IN THE PRESENCE OF XSSR 



Under the spin rotation defined in section V, I would like to underscore the following 
transformation in the equation of motion for the one-particle density matrix 
dH(k, r; T) dHjk, r; T) ^ x dHjk, r; T) 

+ s - (r) ok ■ (A1) 

Because of gauge fields in the rotated basis, the asymmetric component of the density 
matrix acquires an additional term representing the transverse motion of electrons in the 
presence of topological fields. Calculations show 

a„ ^ „ r 9 dH(k,r;T) d x / \ dH(k,r;T) d\ S/1 _ N IK , 

/(k, r; u , T) = r a vn tA -- + + M 8 -^-^1 ® * r). (A2) 

The solution of p s is still given by Eq.ll. 

Taking into account the asymmetric component derived above, I obtain various charge 
and spin currents in both longitudinal and transverse directions. Supercripts l,t stand for 
longitudinal and transverse directions, c and s stand for charge or spin currents; 1, 2 are for 
transport and pumping currents respectively. 

Various charge currents are 
, H/ , rdu d d k _ c dH(k,r;T) d n „ 

(A3, 

Spin currents are 



mi 4vr (2vr) d L ^ v ' dr v 1 de k ' 

J ^' r » -f/l^r^- ' ^ w 8 8 4&>**™m (A4) 
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d = 2, 3 is the dimension of fermi seas. 

Currents J tcl and J z > tsl in these expressions yield contributions to the anomalous Hall 
effect and spin Hall current respectively. Superscript SI has been introduced to specify the 
symmetric part of density matrix which is linear in the external frequency. 



APPENDIX B: CURRENT EXPRESSIONS IN THE PRESENCE OF KSSR 

In the presence of KSSR, the longitudinal spin or charge currents are the same as given 
in the previous subsection. The transverse spin and charge transport (1) or pumping (2) 
currents are 

J * 1(r) = / ^(§5 Tr{n(k) ■ <- ik) ® ?E isr 1 ' f <*»)} 

j rW = / ^(g^Hk) ■ rf£(k) a BH( >£ T) i-Af»'(k, k;„, T)}. (Bl) 



APPENDIX C: ZEEMAN FIELD DEPENDENCE OF THE HALL AND SPIN 
HALL CURRENTS IN THE RASHBA MODEL 

In this section I am going to apply results obtained in Appendix B to briefly discuss the 
Hall current and spin Hall current. The spin hall current in the Rashba Model has been 
studied in various papers and I do not intend to reproduce all results in details. Here I will 
address the effect of Zeeman fields using the approach developed above. 

Define the Hall conductivity and spin Hall conductivity in a usual way, 

3 y = a yx E x , J z y = a z yx E x . (CI) 

After taking into account the symmetry of S^, I find the following expressions for spin Hall 
and Hall conductivity in terms of S^ 2 , 
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^ = ^/?^ 2n2(k)s ^ (k)no(ek) ' 



<V = j ^|2^r Pk [l + / 2 (k)] 1/2 Sf/(k) A no(ek ). (C2) 



2 7 (2tt) 5 

The rest of calculations is straightforward and I derive the following results for the spin 
Hall and Hall conductivity in terms of /(pk) 

lr, I dl 



a. 



yx 



An J rK (l + / 2 ) 2 9p k 



1 r gn B B (pk) dl dn {e k ) 

ayx = J dp * i+p w^^r (C3) 

Finally, I obtain the Zeeman field dependence of the spin Hall and Hall conductivity, 

1 kp 



yx 8n\ 2 s + ky 

1 gfj, B B kp 

ayx = Tn e F Af+fcJ" (C4) 

\ s is defined in section VI B. As mentioned before, I have set % = e = 1 in this article. 

&y X at zero Zeeman field was studied in a few recent works using Kubo formula [12,41]. 
I find that crf„ is a smooth monotonous function of Zeeman fields and decreases as fields 
increase. a yx has a maximum at A s = fc^. Note all results are valid when the adiabaticity 
conditions in section VLB are satisfied. Since the main contribution to the spin Hall current 
is actually from states close to |k| = 0, it is essential to lift the degeneracy at k = 0. So in 
the presence of impurity scattering, the adiabaticity condition implies that a finite Zeeman 
field need to be present and the limit of zero field should be taken with great care. 



39 



REFERENCES 



[1] R.B.Laughlin, Phys.Rev.B 23, 5632(1981). 
[2] D. J. Thouless, Phys.Rev.B 27, 6083(1983). 

[3] Q. Niu and D. Thouless, J. Phys. A 17, 2453(1984); Q. Niu, D.J. Thouless, and Y.S.Wu, 
Phys.Rev.B 31, 3372(1985). 

[4] M. Berry, Pro. R. Soc. Lond. A392, 45(1984). 

[5] B. Simon, Phys. Rev. Lett. 51, 2167(1983). 

[6] D. J. Thouless, M. Kohmoto, P. Nightingale, and M. Den Nijs, Phys. Rev. Lett. 
49,405(1982). 

[7] J.W. Ye, Y. B. Kim, A.J. Millis, B.I.Shraiman, P. Majumdar, and Z. Tesanovic, Phys. 
Rev. Lett. 83, 3737(1999); Y. B. Kim, P. Majumdar, A.J.Millis, and B.I.Shraiman, 
cond-mat/9803350. 

[8] P. Matl, N. P. Ong, Y. F. Yan, Y. Q. Li, D. Studebaker,T. Baum, and G. Doubinina, 
Phys. Rev. B57,10248(1998). 

[9] G. Sundaram and Q. Niu, Phys. Rev. B 59, R12755(1999). 

[10] T.Jungwirth, Q. Niu and A. H. MacDonald, Phys. Rev. Lett. 88, 207208(2002). 

[11] R. Karplus and J. M . Luttinger, Phys. Rev. 95, 1154 (1954). 

[12] J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, and A. H. MacDonald, cond- 
mat/0307663; D. Culcer, J. Sinova, N. Sinitsyn, T. Jungwirth, A. H. MacDonald, and 
Q. Niu, cond-mat/0309475; N. Sinitsyn, E. M. Hankiewicz, W. Teizer, and J. Sinova, 
cond-mat/0310315. 

[13] M. Onoda and N. Nagaosa, J. Phys. Soc. Jpn.71, 19(2002); G. Tatara and H. Kawamura, 
J. Phys. Soc. Jpn. 71, 2613(2002). 

40 



[14] S.Murakami, N. Nagaosa, and S. C. Zhang, Science 301,1348(2003). 

[15] S.Murakami, N. Nagaosa, and S. C. Zhang, cond-mat/0310005. 

[16] B. A. Bernevig, J. P. Hu, E. Mukamel, and S. C. Zhang, cond-mat/0311024. 

[17] F.Zhou, B. Spivak, and B. Altshuler, Phys. Rev. Lett. 82, 608(1999). 

[18] F.Zhou, cond-mat/0202321; see also F. Zhou, Phys. Rev. B 65, 220514(2002). 

[19] J. M. Luttinger, Phys. Rev. Rev. 102, 1030(1956). 

[20] G. Dresselhaus, Phys. Rev. 100,580(1955). 

[21] R. I. Rashba, Soviet Physics Solid State 2, 1109(1960). 

[22] E.I. Rashba and V. I. Sheka, Sov. Phys. Solid State 3, 1257(1961). 

[23] Y. A. Bychkov and E. I. Rashba, Sov. Phys. JETP Lett. 39,78 (1984). 

[24] E. I. Rashba, cond-mat/03090441. 

[25] The notion of meron was first introduced in D. J. Gross, Nucl. Phys. B 132, 439(1978). 
I want to thank Ian Affleck for bringing this reference to my attention. 

[26] The mass generation by merons in spin chains was discussed in I. Affleck, Phys. Rev. 
Lett. 56, 408(1986). See also related discussions on spin gap in F. D. Haldane, Phys. 
Rev. Lett. 50, 1153 (1983). 

[27] B. Spivak, F. Zhou, and M. T. Beal Monod, Phys. Rev. B51, 13226(1995). 

[28] This is a stringent sufficient condition. However if the main contribution to currents is 
from electrons with momentum kp away from the k = point, the condition can be 
relaxed and is always satisfied. 

[29] When the contributions to currents are mainly from the Fermi surface as in the toy model 
discussed in section VI. A, the adiabaticity condition can be easily satisfied because of 



41 



the finite splitting between two spin bands due to spin-orbit coupling. However, in the 
2D models studied here contributions might actually come from states in the vicinity 
of k = 0; to establish the adiabaticity condition it is essential to introduce a finite 
Zeeman field so to split the degeneracy at k = as shown in Eq.100. Then the Born 
approximation is still applicable when the impurity potentials are weak compared with 
the finite Zeeman field splitting. 

[30] P. Brouwer, Phys. Rev. B58, 10135(1998). 

[31] F.Zhou, Int. Jour. Mod. Phys.B, 117(2001). 

[32] M. Switkes, C.Marcus, K. Campman, and A. Gossard, Science 283, 1907(1999). 
[33] T. A. Shutenko, I. Aleiner, and B. I. Altshuler, Phys. Rev.B 61, 10366(2000). 
[34] L. S. Levitov, cond-mat/0103617. 

[35] P. Sharma, and C. Chamon, Phys. Rev. Lett. 87,096401(2001). 

[36] R. Mucciolo, C. Chamon, and C. Marcus, Phys. Rev. Lett. 89, 146802(2002). 

[37] J. Wu, B. Wang, and J. Wang, Phys. Rev.B66, 205327 (2002). 

[38] W. Zheng, J.L. Wu, B. G. Wang, J. Wang, Q. F. Sun, and H. Guo, Phys. Rev. B68, 
113306(2003). 

[39] S. K. Watson. R. M. Potok, C. Marcus, and V. Umansky, cond-mat/0302492. 
[40] P. Sharma and P. Brouwer, Phys. Rev. Lett.91, 166801(2003). 
[41] J. Schliemann and D. Loss, cond-mat/0310108. 



42 



FIGURES 




FIG. 1. Comparison of the usual Stern-Gerlach beam splitting and the topological 
Stern- Gerlach splitting. When the Zeeman fields are applied along the z-direction as shown in 
Fig. (a), in the conventional Stern-Gerlach splitting S z is a good quantum number and no spin 
rotation occurs; the Zeeman field gradient drives spin- up and spin-down particles apart along the 
y-direction. However TS GS is always accompanied by spin rotation as shown in Fig.(b). Spins are 
represented by short arrows in this Figure and other figures. 
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FIG. 2. Plus and minus Fermi seas for the toy Hamiltonian in Eq.65, the Rashba Hamiltonian 
and Dresselhaus Hamiltonian. Fig.(a),(b) are for the toy Hamiltonian. In the plus Fermi sea each 
spin points at the direction of its momentum (shown in Fig. (a)); in the minus Fermi sea shown 
in Fig.(b), the spin of an electron points at the opposite direction of its momentum. In Fig. (c) 
and (d), we show the corresponding Fermi seas for the 2D Rashba Hamiltonian. In k x — k y plane, 
electron spins in the plus fermi seas form a meron with a half skyrmion charge when a Zeeman 
field is applied along the z-axis. At large k limit in the plus Fermi seas, spins point at e x O(k) 
direction while in the minus Fermi seas spins point at — e x fi(k); they both represent vortices 
with one unit positive vorticity. Here fi(k) is a unit vector along k. At the center of fermi seas, 
spins are along the ±z directions (not shown here). In Fig.(e) and (f), we show spin rotation in 
fermi seas of the Dresselhaus model; electron spins in the plus fermi seas form a meron with a half 
negative skyrmion charge. At large momentum, electron spins in both fermi seas form vortices 
with one unit negative vorticity. 




FIG. 3. Topological field T as a function of /?k in the 2D Rashba model (schematic). X s which 
defines the size of meron is chosen to be smaller than the fermi momentum fcp 
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FIG. 4. Responses of Fermi seas and individual electrons to external pumping fields applied 
along the x-direction. Fig. (a) is for polarized electrons in the presence of orbital magnetic fields 
directing along the z-direction; spin-up and spin-down Fermi seas experience identical displacement 
d along the longitudinal direction x and transverse direction y which are shown in two upper left 
insets of Fermi seas. Fig. (b) is for electrons with spin rotation in the X-space; in this case 
spin-up and spin-down Fermi seas have opposite displacement d along the transverse direction y as 
shown in two insets for two Fermi seas. In Fig.(c) we illustrate differences in electrons' responses 
in the plus and minus Fermi seas. We want to emphasize that Fermi seas don't have collective 
displacement along a transverse direction y or d^ = 0. In insets we show the external-field-induced 
drift or group velocity of individual electrons v y in different regions of fermi seas. The transverse 
drift of an electron in the upper part of a Fermi sea with momentum k is in an opposite direction of 
the drift of the electron with momentum — k in the lower part, which results in self- twists of Fermi 
seas. In upper(lower) insets, we show the drift of electrons in the north (south) poles of plus and 
minus fermi seas respectively; the big arrows across Fermi seas indicate the direction of two distinct 
twists of plus and minus Fermi seas. In Fig.(d), we illustrate responses of the plus and minus fermi 
seas in the 2D Rashba model; in this case, again two fermi seas have zero displacement in the 
transverse y-directions (d y = 0). However, the spin-plus and spin- minus electrons in two fermi 
seas acquire a field-induced dispersive group velocities in the opposite y-direction as shown in two 
insets in Fig.(d); the field- induced transveral group velocity decreases rapidly as the momentum 
increases. See section VII for detailed discussions. 
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